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Abstract 

A multi-species generalization of the Asymmetric Simple Exclusion Process (ASEP) has been considered 
in the presence of a single impurity on a ring. The model describes particles hopping in one direction with 
stochastic dynamics and hard core exclusion condition. The ordinary particles hop forward with their char- 
acteristic hopping rates and fast particles can overtake slow ones with a relative rate. The impurity, which 
is the slowest particle in the ensemble of particles on the ring, hops in the same direction of the ordinary 
particles with its intrinsic hopping rate and can be overtaken by ordinary particles with a rate which is not 
necessarily a relative rate. We will show that the phase diagram of the model can be obtained exactly. It 
turns out that the phase structure of the model depends on the density distribution function of the ordinary 
particles on the ring so that it can have either four phases or only one. The mean speed of impurity and 
also the total current of the ordinary particles are explicitly calculated in each phase. Using Monte Carlo 
simulation, the density profile of the ordinary particles is also obtained. The simulation data confirm all of 
the analytical calculations. 
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1 Introduction 

One of the simplest models of driven diffusive systems is the Asymmetric Simple Exclusion Process (ASEP) 
§ § |, §. The ASEP consists of particles on a one-dimensional lattice % attempting to jump to an adjacent 
site, choosing the site to their right with rate p, and to their left with rate q, with the attempt succeeding if the 
target site is not already occupied. This model can describe the Burgers equation on a microscopic scale 
However the ASEP is not only interesting because of its relation to the Burgers equation but also it has many 
applications in other fields such as growth models, traffic and kinetics of portein synthesis [0, ^, [|. During 
last decade the ASEP has been studied extensively both with open boundaries and periodic boundary condition. 
It is known that this simple model shows a variety of interesting phenomena such as boundary-induced phase 
transitions and shock profiles Jh], [ll], [lj, [b|] . 

In the open boundary case the particles are injected and also extracted from both ends of the chain so that 
the density of particles on the chain might vary. In this case the density of the particles on the chain can be 
controlled by some boundary parameters, which in fact control the input and output rates. Depending on the 
value of these parameters, the phase diagram of the model can possesses three different phases: a high-density 
phase (HD), a low-density phase (LD) and a maximal-current phase (MC). It is known that the density profile of 
the particles has a shock structure ^ on the coexistance line between LD and HD phase where a first order phase 
transtion takes place. A second order phase transition also exists from LD and HD phases to the MC phase 
|l5| , |l6| . This phase structure seems to be generic for all those models in which the fundamental digram (the 
current of particles as a function of the mean density of them) has a single local maximum. The ASEP with 
periodic boundary condition, where the total number of particles on the chain is a conserved quantity, in the 
presence of a second class particle has also been studied in two different cases |l7], [[ll |l9|, |2(| ^l) . In the first 
case the second class particle (which is called the impurity) hops only to the right. The ordinary particles hop 
both to the left and right and can also overtake the impurity from the left. The phase diagram of the model in 
comparison to the open boundary problem is rather complicated. Again the system presents a shock in one of 
the phases. In the second case the impurity hops only to the left. The ordinary particles and the impurity may 
exchange their positions with a finite rate. This model has only been studied for q = 0, and the case q ^ is 

2 The "shock structure" means a sharp discountinuty between a region of high density of particles and a region of low density. 
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still an open probelm. It has been shown that the phase diagram of the model for q = has only two different 
phases: A shock phase and a low-density phase In the absence of the second class particle the stationary 
state is well known: all configurations are equally likely. 

Recently, a multi-species generalization of the ASEP which consists of p different species of particles hopping in 
one direction has been introduced [2^| . Since the number of species is quite arbitrary, one can take the infinite 
species limit. In this limit the hopping rate of each species of particles v is taken from a continuous distribution 
<j(v). In this model the overtaking between particles is allowed so that the fast particles can overtake the slow 
ones. Particles of different kinds can also be inserted and extracted from the system. Using a Matrix Product 
formalism (first introduced in [p^[) the phase diagram of this model in the stationary state has been obtained 
p4[ . It turns out that this multi-species generalization has a richer phase structure (in comparison to the usual 
ASEP) depending on the distribution function a(v) so that it can exist either in a three-phase or a two-phase 
regime. The phase diagram of the model in the three-phase regime is nearly similar to the phase diagram of 
the ASEP with open boundaries. However, in the two-phase regime the high density region disappears. 
In this paper we study the multi-species ASEP in the presence of an impurity on a ring. We will consider the 
case in which the impurity moves in the same direction of the ordinary particles and can be overtaken by them. 
In the periodic boundary condition case, instead of the distribution function o~(v) which gives the probability of 
finding a particle with hopping rate v before entering the system, we introduce p(v) which gives the density of 
the ordinary particles with hopping rate v on the ring. This quantity can also give the probability of finding an 
ordinary particle with intrinsic hopping rate v on the ring at any time. We will show that the phase structure 
of this model depends on the density distribution function of normal particles p(v) on the ring so that it can 
be in two different regimes. In the first regime the phase diagram is rather similar to the phase diagram of the 
one-species ASEP with an impurity on a ring (the first case) |Q and it has four different phases. There exists 
a region in which the total density profile of the ordinary particles has a shock structure. In the one-species 
limit the speed of impurity in the shock phase is independent of the total density of the ordinary particles on 
the ring. However, for our multi-species model the speed of impurity in this phase, which is in fact equal to the 
speed of the shock, depends on p(v). In the second regime, we may see two different behaviours. If the speed of 
impurity is smaller than the hopping rate of the slowest ordinary particle on the ring, then the phase diagram 
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will contain only one phase. However, if we allow the impurity have a hopping rate greater than the hopping 
rate of the slowest ordinary particle, then a new phase will appear so that the phase diagram in this regime will 
contain two different phases. In order to check the analytical results we have done the Monte Carlo simulations. 
It turns out that the simulation data confirms all of the analytical results. Our analytical approach does not let 
us calculate the density profile of the ordinary particles on the ring; however, we will try to obtain the density 
profile of the ordinary particles on the ring using the Monte Carlo simulations. The multi-species ASEP has 
also been studied on a ring in the absence of the impurity p3| . It is known that the stationary state of this 
model is uncorrelated. We will show that on an specific line in the parameters space the impurity behaves like 
an ordinary particle so that on this line the same results of [^3| can be derived. 

This paper is organized as follows. In section 2, we will define the model and obtain its generating function for 
the grand canonical partition function on the ring. Analyzing the singularities of this generating function lets us 
find the phase structure and also calculate the mean value of two important physical quantities, i.e. the mean 
speed of the impurity and the total current of the ordinary particles in the stationary state. In section 3, we will 
study two different examples to check the validity of our calculations: First we will take the one-species limit 
and show that all of our results converge on those obtained in fl§| . In the second example we will introduce a 
family of density distribution functions, and using Monte Carlo simulation show that all of the analytical results 
obtained for this family of distribution are correct. In the last section we will present the conclusions. 



2 Exact Stationary State Using the Matrix Product Formalism 
2.1 Definition of the Model 

Our model consists of p species of first class particles (we refer to them as the ordinary particles) and a single 
second class particle (we refer to it as the impurity) on a chain of length L + 1 with periodic boundary condition. 
We assume that the number of particles on the ring is a conserved quantity. Each species of particles has its 
own intrinsic hopping rate. Fast ordinary particles overtake the slow ones with their relative hopping rate. The 
impurity has the smallest hopping rate in the ensemble of the particles on the ring and thus, can be overtaken 
by ordinary particles of any kind with a rate which is not necessarily a relative rate. In this paper, following ]24j]. 
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we take the infinite species limit, since the number of species is an arbitrary number. In this limit the hopping 
rate of the ordinary particles v is a continuous parameter and the density of particles in each species has also 
been taken from a continuous density distribution p(v). The function p{v) is defined in the interval [v%, +oo( in 
which vi is the smallest hopping rate in the ensemble of the ordinary particles. Denoting an ordinary particle 
with the intrinsic hopping rate v by D(y), the impurity by A and an empty site by E, during the infinitesimal 
time step dt any bound + 1) (with 1 < i < L + 1) evolves as follows: 

D(v)E — > ED(v) with the rate v 

D(v)D(v') — ► D(v')D(v) with the rate v - v' , v > v' (1) 

D(v)A — > AD(v) with the rate v + 0-1 

AE — > EA with the rate a 

in which a and /3 are positive parameters. We specify each ordinary particle by its intrinsic hopping rate 
hereafter. According to the above process, particles hop to their rightmost site with their intrinsic hopping rate, 
provided that it is empty. If an ordinary particle of kind v encounters a site occupied by an ordinary partice of 
kind v' , they might exchange their positions with rate v — v' provided that v > v' . An ordinary particle of kind 
v can also overtake the impurity with rate v + (3 — 1. For this rate to be positive, we require that /3 > 1 — V\. 
It can be seen that, as far as a + (3 ^ 1 the impurity behaves quite differently from the ordinary particles. For 
a + j3 — 1, since a < V\, the impurity behaves like an ordinary particle with the smallest hopping rate. On this 
line, all the ordinary particles overtake the impurity with a relative rate. It is known that the steady state of the 
system in this case is completely uncorrelated and all the configurations have the same probability of occuring 
[ p3[ . Also, the special case p(v) = p5(l — v) has already been studied in Jig ]. In this case all of the ordinary 
particles have the same hopping rate (one-species limit) so that no overtaking between them takes place. In 
order to study the stationary state properties of this model, we will introduce the necessary mathematical tools 
in the following section. 
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2.2 Algebraic Preliminaries 

Among different methods the Matrix Product formalism has proved to be one of the most useful for obtaining 
the stationary state properties of reaction-diffusion models. According to this formalism the stationary proba- 
bility P({C}) of any configuration {C} can be written as a trace of the product of non-commuting operators 
[ p5| . In our model, for the sake of simplicity, we will work in the relative frame of the impurity i.e. using the 
translational invariance of the system we always keep the impurity at site L + 1. We will also use a grand 
canonical ensemble, in which the number of ordinary particles in each species fluctuates arround a mean value. 
However, the density of each species should be fixed by its fugacity. By using the MPA the stationary probability 
of any configuration, {C} can be written as follows: 

P({C}) = ±Tr(f[x i ) (2) 

i— 1 

in which = E (1 < i < L), if site i is empty, and = z(v)D(v) (1 < i < L) if it is occupied by an ordinary 
particle of kind v. Since the site L + 1 is always occupied by the impurity, we take X^ + \ — A. The quantity 
z(v) is the fugacity of the ordinary particles of kind v which should be fixed using: 

z ( v^) S 

p(v) = UnfiL^oo— — — —lnZ L . (3) 
L oz(v) 

We have not introduced the fugacity of holes since the total number of available sites L + 1 is fixed. The 
total density of the ordinary particles is also specified by J p(v)dv = p. The normalization factor Zl in the 
denominator of (||) plays the role analogous with the grand canonical partition function in equilibrium statistical 
mechanics, and can be calculated using the fact that J^c -^({C}) = 1- Thus one finds: 

Z L = Tr{C L A) , C = E+ f z(v)D(v)dv. (4) 

As we will see shortly, the mean value of the physical quantities such as the speed of impurity and the current 
of each species of ordinary particles can be written in terms of Z^. The operators E, D(v)'s and A then satisfy 

3 This formalism is also known as the Matrix Product Ansatz (MPA) in related literature. 
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the following relations: 

D(v)E = ±D(v) + E 

D(v')D(v) = ^r(vD{v') - v'D(v)), v'>v 

(5) 

D(v)A=^ =I A 
AE = -A 

a 

Following j26) we write A as |V)(W| which simplifies the above algebra to the one in [^4| for the study of the 
multi-species ASEP with open boundaries. 

In the following we will show how one can calculate the mean value of physical quantities in the stationary state 
using the MPA. For example, the speed of impurity in the reference frame of the lattice can be written as: 

Tr{C L ' x AE) f. a 1 .Tr(C I - l z(v)D(v)A) , 
V = a Tr(C L A) -] V } V + P - l) Tr(C L A) ^ 

Using the algebra of (H) one can easily simplify the above expression to find: 



V = (l- I vz{v)dv) T j}- G .^ T 1 f > = (1 - / vz(v)dv)^^. (6) 



One can also find a simple expression for the total current of the ordinary particles on the ring Jtotai in the 
stationary state as follows: We specify each configuration of the system by an L + 1-tuple (n, • ■ • , tl+i) where 
Tj = v if site i is occupied with an ordinary particle of kind v , Tj = a if it is occupied by the impurity, and Tj = e 
if it is empty. Now one can write the current of the ordinary particles of kind v, J(v), in the stationary state 
as follows: 

J(v) = vP(r k = v,r k+ i = e) + j v , <v {v - v')P{T k = v,T k+ i = v')dv' 
- L>>v( v ' - v ) P ( T k = v '> T k+i = v)dv' + {v + f3- l)P(n = v, T k+1 = a) 
where P{r k = ■ ■ ■ , T k +\ = • • •) is the probability of finding the system in configuration {C} provided that the 
sites k and k + 1 are occupied by particles or holes. The above relation in turn can be written in terms of 
the conditional probability P(r k = ■ ■ ■ ,T k +i = ■ ■ ■ \ti = ■ ■ ■) given by (||) which has the same definition as 
P(r k = ■ ■ ■ , Tfc + i = • • •) provided that the impurity is at site /: 

J(v) = vYsiti 1 P( T k = v,r k+1 = e\n = a)P(ri = a) + 
T,iti 1 L> <v ( v - v')P{T k = v,r k+1 = v'\n = a)P(n = a)dv' 
~ Ya=i L>>v( v ' ~~ v )P( T k = v', r k+1 = v\n = a)P(n = a)dv' + (v +(3 - l)P(r fe = v, r k+1 = a). 



where the prime over a sum means that the terms I = k, k + 1 should be excluded. Since the system in 
the stationary state is translationally invariant, the impurity can be on any of the L + 1 sites with the same 
probability. Now using the algebra (JsJ) we obtain: 

1 vz(v) TrUJ^A) 
P(r k = v,r k+1 =a) = P(r k = v\r k+1 = a)P{r k+1 = a) = ——^ __ 

which after calculations, results in the following expression for the current of the ordinary particles of kind v: 

J( v ) - v ~( r ] L TrjC^A) _j_ V L-1 TrjC"- 1 z(v)D(v)C^- k ^ A) 
J(U) — UZ(U) L+1 Tr ( C L A) r L+l 2—ik=\ Tr(C L A) 

vz \ v ) L+l 2^ k =l J Vl Tr{C L A) UU + 



f z(v)z{v')(v-v') Tr(C k - 1 D(v)D(v')C L - k - 1 A) j ,_ 

l^k=l Jv'<v L+l Tr(C L A) aV 

r z(v)z(v')(v'-v) Tr(C fc ~ 1 D(v')D(v)C L ^ k ' 1 A) , / 

2^fc=l Jv>>v L+l Tr{C L A) av ' 

Using the definition of the total density of the ordinary particles of kind v (||): 

p{v) = I ^ t^Wa) 

and also the algebra of (|^) one can find the following simple expression for J(v): 
m(L'<,-L'>J(^)Mvl-z(v')v'p(v))dv']^^ 



In the thermodynamic limit L — > oo we find: 

J(v) = [p(v) + (1 - p)z(v)v + ([ - f )(z{v)vp(v') - ziv'ypiv))^'} 21 - 1 



Z L 

In order to calculate the total current to the ordinary particles J total , one should integrate the above expression 
over all species: 

Jtotal = J Vl J{v)dv = ^ 

pV + (J vi vz(v)dv)^ =p V +(^- V). 

As can be seen both Jtotal and V are functions of Jr 1 and also independent of the position on the ring as 
one should expect in the steady state. In principle, one can calculate these quantities using the non-trivial 
(non-commuting) representations of the algebra (|^) . It has been shown that all non-commuting representations 
of (||) are infinite dimensional |22| . However, the author of has shown that a coherent basis exists which 
using it one can compute the generating function for the canonical partition function of this model with open 
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Figure 1: Plot of S(z; it) (Left) as a function of u in the interval u £ [0, ui] and also B(z; s) (Right) the inverse function 
of S(z; u) in the interval s £ [0, s m ] as a function of s. 

boundaries p4j. Finding the singularities of this generating function lets us obtain not only the exact phase 
diagram but also the mean values of the physical quantities. The generating function for the grand canonical 
partition function of our model which is defined as: 



/(z;,s;a,/3) = £ s L Tr(C L A) = £ s L (W\C L \V) = £ s L Z L . 



(8) 



L=0 



L=0 



L=0 



can also be calculated exactly using the same coherent basis. In fact (g) has already been calculated to obtain 
the mean density of the ordinary particles in the middle of chain in the open boundary problem (see equation 
(64) of p3|) and is given by the following exact expression: 



in which: 



/(z;s;a,/3) = 



S{z;u) 



uB(z;s) _ q(l-j3) 
a-B(z;s) a+P—1 

(l-^)[^f^j-ll 



, H6[0,«i]. 



(9) 



(10) 



I + f £14^ 

In (|^) the function B(z; s) is the inverse of S(z; u) in the interval u 6 [0, u m ] which tends to zero as s — > and 
also has the following properties: 

( S[z;B(s)] = s, 0<s<s m 

{ (11) 
[ B[z; S(u)] — u, < u < u m . 

In figure 1 we have plotted S(z; u) and also B(z; s) as a function of u and s respectively. One can easily see 
from ( |Io|) that for u £ [0, Vi], -j^zS^ 1 is always positive so S~ x is a concave function in that region. On the 
other hand since z(v) is a positive real number, ( |ic| ) is a positive function of u for u G [0, V\\ so it can only have 
one maximum u m in this interval. The existence of the local maximum depends on the sign of 



du 



S 



u—vi V\ 



vz{v) 
V1 (v-vif 



-dv. 



(12) 



Because of the relation (|3j) , the above quantity in turn is a function of the density distribution function of the 
ordinary particles on the ring p(v). The generating function ^) is restricted to the region a + (3 > 1 so in order 
to have access to the whole parameters space we will have to introduce an analytical continuation. 
In the next section we will obtain both the phase diagram and also the values of V and Jtotai by analyzing the 
singularities of ([))). 

2.3 Exact Results for the Physical Quantities and the Phase Diagram 

As we mentioned, the phase structure of the model and also the mean value of the physical quantities are 
determined by the singularities of (JsJ) . The number of singularities of (^) depends on whether the local maximum 
exists or not. If ([l2]) is negative, then the local maximum exsists, and (|J) has four singularities depending on 
the values of a and 0: 

s a = S(z; a) for a < u m , S(z; a) < S(z; 1 — 0) (13) 
s$ = S(z; 1-/3) For < 1 - u m , S(z; l-0)< S(z; a) (14) 
s m = S(z; u m ) for a > u m , > 1 — u m (15) 
so = s a — s/3 for a < u m , < 1 - u m . (16) 

The first (second) singularity arises from the disappearance of a — B(z;s a ) ( s(*-i-p) ~ ^) ^ n numera tor 
(denominator) of (^). The third singularity occurs at s m where S(z;u) has a local maximum at u m , and so B 
becomes singular. In order to calculate u m in terms of the density distribution function p(v), we note that the 
convergence radius of the formal series (|^) can be written as: 

R(z)= lim ((W\C L \V))~* (17) 

so the mean density of the ordinary particles of kind v in each phase (^|) can be written in terms of the 
convergence radius R(z): 

S 1 

p(v) = z{v)— T ln— . (18) 
oz(v) R[z) 

On the other hand, the convergence radius of (^|) is the absolute value of its nearest singularity to the origin; 
therefore, from (jlfl), ( |l5| ) and (|l^) we arrive at the following expression which gives u m only as a function of 
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p{v): 

vp{v) 



■dv = l. (19) 

The last singularity takes place whenever s a = sp. Using ( |Io| ) one finds the following condition on z(v), a and 
f3 in this phase: 

vz(v)dv 1 (2Q) 



l Vl (v-a)(v + P-l) a(l-0) 
The above condition, which is the consequence of using the grand canonical ensemble, means that only two of 

the mentioned quantities can be chosen freely. We consider two different cases: First, we assume that (3 is a 

dependent variable. By using ( |l0| ) and ( Ug ) we obtain: 

p(y) = a(l-p)p&. (21) 



Substituting @ in @ gives: 

vp(v)dv 



1. (22) 



Ul v + - 1 

Now comparing ( |l9| ) and (p3) this specifies one of the boundaries of this phase i.e. u m = 1-/3. Similarly, taking 
a as dependent variable, one finds u m = a which gives the other boundary of this phase. In the open boundary 
problem sq determines a line in the parameter space. However, in our model it specifies a region in which, as 
simulation results show, a shock profile evolves in the density profile of the ordinary particles on the ring. In 
order to determine the boundaries of different phases one should note that using ( [Tol) and (|l8| ) the condition 
s a < sp in dig] ) can be written as (3 > 1 — u m . Similarly one can see that the condition sp < s a in (|l4| ) can be 
written as a > 1 — u m . Unforunately, we have not been able to find a simple expression for the total density 
profile of the ordinary particles on the ring defined as: 

f . .j f TrfC 1 ^ 1 z(v)D(v)C L ~ i A) , 

Ptotal(i) = J P (v,i)dv = J Tr{ CL A ) dv ^ 

where p(v, i) is the density of the ordinary particles of kind v at site i. However, we can find the exact values 
of Ptotai{i) exactly in front of the impurity ptotal(l), and also behind it p total (L). For p tota i(l) we obtain: 

f Tr{C L ^Az{v)D{v)) 1 Tr{C L ^A) 

PtoUl) = dv = 1 - - Tr(cM) (24) 

in which we have used the algebra (^|). Also for ptotai(L) we obtain the following expression using the algebra 



PMiL) . , za^fe^a* . ( r j^lj^^mi. (25 , 



Tr(C L A) y J Vl v + 13-l' Tr(C L A) 

10 



p 



1-V, 



Figure 2: Phase diagram of the model where -^S 



< 0. The short lines show the behaviour of the total density 



profile of the ordinary particles obtained from simulation results. The dotted line is a + j3 = 1. 



Using @ and (|T|), one finds p to t a i(l) = p and also p to tai{L)\ L ^ oo = p- (a + (3 - 1) ^^ dv in the phase 
( |l3|) . Note that while the total density of the ordinary particles is constant directly in front of the impurity, it 
might be greater (if a+(3 < 1), smaller (if a+(3 > 1) or equal (if a+(3 = 1) to p immediately behind it. The same 
phenomenon might take place with a slight difference for the total density of the ordinary particles immediately 
in front of the impurity in the phase jl^ ) where we obtain ptotaliX) = P+ ( a+ ^~ c ^( 1 ~ p ) and ptotal{L)\ L ^ — p- 
In this case ptotai{^) can be greater (if a + (3 > 1), smaller (if a + (3 < 1) or equall (if a + j3 = 1) to p. In order 
to see the behaviours of the total density profile of the ordinary particles on the ring, we have performed the 
Monte Carlo simulation. We will discuss the obtained results in the next section. 



In figure 2 we have plotted the phase diagram of our model for -^S 



< 0. The total density profile of the 



ordinary particles (23) obtained from the Monte Carlo simulation is shown as a curved line in each phase. One 



should note that ptotai(i) has two different behaviours in the phase (13). The same phenomenon can be seen in 
the phase (|l4|) (see figure 2). As we have already mentioned, the line a + [3 = 1 is where the impurity behaves 
like the slowest ordinary particle. Simulation results show that the total density profile of the ordinary particles 
on this line is flat. These results are in agreement with the picture that one can obtain from our analytical 
results for the values of the total density of the ordinary particles immediately behind as well as in front of 
the impurity. In the region specified by ([l6]) the single impurity provokes a macroscopic shock in the system. 
The density profile of the ordinary particles in the shock phase, obtained from the Monte Carlo simulation, can 
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1-V, 



Figure 3: Phase diagram of the model where 4-S 



> 0. The short curved lines show the bahaviour of the total 



density profile of the ordinary particles obtained from simulation results. The dotted line is a + /3 = 1. 

be seen in figure 4. Nevertheless, when (|l^) becomes positive, the function ([h]) is a monotonically increasing 
function of u and the formal series (0) has only one singularity at s a : 



s a = S(z;a) for a < v\. 



(26) 



This singularity comes from the disappearance of 1 — aB(z; s a ) in the numerator of (|^). The simulation data 
predicts that in this case, the total density profile of the ordinary particles (^3|) has two different behaviours in 
the parameters space. Again, on the line a + (3 — 1 the total density profile in the steady state is flat. In figure 
3 we have plotted the phase diagram of our model for 4-S > 0. 

u— vi 

The physical quantities V and J total , which are given by (^J) and (0) respectively, can be written in terms of the 
convergence radius of (0) : 



V = (1 - J vi vz(v)dv)R(z) 

Jtotal =pV+ (R(z) - V) 

By using (IS) and (E7]) one obtains the following expressions for V and Jtotal in each region of (n3|)-(nq): 



(27) 



V = a - f Vi vp{v)dv , Jtotal = (1 - p) J Vl vp(v)dv for a < u m , (3 > 1 - u r , 

V = 1 - /3 - J vi vp(v)dv , J to tai = (1 - p) J V1 vp(v)dv for a > u m , /9 < 1 - u r , 

V = u m ~ J vi vp(v)dv , Jtotal = (1 - p) J Vl vp(v)dv for a > u m , fi > 1 - u r , 

V = (a- P + l- u m ) - f Vi vp(v)dv , Jtotal = ( ° „~ - V)(l - p) for a < u m , p < 1 - u r , 
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(28) 



As can be seen from (Pq ) both the speed of impurity and total current of the ordinary particles are continuous 
functions of a and /3. We should note that using the grand canonical ensemble procedure, one cannot find the 
exact expressions of V and Jtotai i n the shock phase except on the boundaries of this phase. This phenomenon 
has already been observed in [|l9[ |o) ; nevertheless, the Monte Carlo simulation data in this phase predict that 
for a fixed value of a {(3) these quantities are linear functions of (3 (a). Keeping in mind that these quantities 
change continuously over the boundaries, it is not difficult to conjecture the expressions given in the shock phase 

however, Monte Carlo results confirm the correctness of these expressions. 
Using (^6|), (|l8|) and (p7|), we also obtain the following expressions in the second regime: 

V = a- I vp(v)dv , J to tai = (1 — p) \ vp(v)dv for a < v\ (29) 

J vi J vi 

Remark: 

If we had released the condition a < v\ (the fact that the impurity should be the slowest particle on the ring), 
there would have found a change in the phase diagram of our model in (p6|). It can be shown that in this case 
an extra phase exists for a > v x which is produced by the singularity of (|l0|) at v\ 0: 

s Vl = S(z;vi) for a>v\. (30) 

The speed of the impurity as well as the total current of the ordinary particles can also be obtained in this 
phase: 

V = Vi- / vp(v)dv , Jtotai = (1 - p) \ vp(v)dv for a > V\. (31) 

Jvi Jvi 

3 Special Cases and Simulation Results 

In this section we will check the validity of our analytical results obtained in the previous section by considering 
two different examples. As the first example we consider p(v) = pS(v — 1). In this case all the ordinary particles 
have the same hopping rates equal to the unity, and p is their total density. This model has already been 
studied in [fl8[ as a simple exclusion model for the study of shocks in one-dimension. The function (|Io| ) has a 
very simple form in this case: 

S(z;u) = -^4-. (32) 
1 + (z — l)u 

4 In fact (^) has a branch cut on a segment of the real axis begining from u = vi to +oo. 
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One can easily see that ( |l2|) is equal to — oo and also ( |l9| ) gives u rn = 1 — p; therefore, the mean values of the 
speed of impurity and total current of particles should be obtained from (GS 



V = a - p , Jtotai = p(l - p) for a<l-p,/3>p 

V = 1 - /3 - p , Jtotai = p(l - p) for a > 1 - p, p < p 

(33) 

V = 1— 2 — p , Jtotai = p(l - p) for a > 1 - p, p > p 

V = a - P , Jtotai = p(a - P) + P(l - a) for a<l-p,P<p 

The results ( p3| ) agree with those obtained in Jl8| ]; however, the analytical calculations in jl8j show that the 
density profile of the ordinary particles in the first and the second phase of ( |33| ) has an exponential behaviour 
with two different characteristic length scales for a + P > 1 and a + P < 1. It is also known that the fourth 
region in ( ]3^ ) is where the density profile of the ordinary particles has a shock structure. 
As the second example, we consider the following density distribution function: 

= {1 p{ Xr^ [v ' vl)n • n -° (34) 

in which n is not necessarily an integer. This kind of distribution has also been considered in related contexts 
|27| |4|. Using |[|) and @ the condition Q can be written as: 



du 



oc ^ - 4f^. (35) 
Vi n{l — vi) 



For n > 1 _ P pl vl the sign of (|35|) becomes positive and the phase diagram should be obtained from ( |26| ) . Taking 
n = 1 and p = 0.6 we will obtain the total density profile of the ordinary particles in each phase using the Monte 
Carlo simulation for v\ — 0.6 and v\ = 0.2. The values of u m for v\ > 0.25 should be obtained from ([l9|). For 
p = 0.6 and v\ = 0.6, using ( p"9| ) and (|34|), we obtain w m = 0.338817. In figures 4 and 5 we have plotted the 
total density profile of the ordinary particles ptotai{i) obtained from the Monte Carlo simulation (as seen by the 
impurity) on a chain of length L = 500 for v\ = 0.6 at a = 0.2 and a — 0.5 respectively. Let us first explain the 
case a = 0.2. Depending on the value of P, three different phases can be seen: a shock phase, where the total 
density of the ordinary particles has a shock structure, and two other phases in which the total density of the 
ordinary particles in the middle of the chain is equal to p. It is also seen that the impurity produces a small 
disturbance just behind itself in these two phases. We have plotted ptotaiif) for oi — 0.5 in figure 5. It can be 
seen that the total density profile of the ordinary particles has again three different behaviours depending on the 
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Figure 4: Plot of the total density profile of the ordinary particles on a ring of length L=500 for a— 0.2, v\=0.6, p = 0.6 
and different values of /3. 
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Figure 5: Plot of the total density profile of the ordinary particles on a ring of length L=500 for a — 0.5, Vi = 0.6, 
p = 0.6 and different values of /3. 
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Figure 6: Plot of the speed of the impurity as a function of j3 obtained from Monte Carlo simulation on a ring of length 
L—500 for vi = 0.6, a — 0.2, a = 0.5 and p — 0.6. The phase transition takes place at /3 C ~ 0.66 
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Figure 7: Plot of the total current of the ordinary particles as a function of obtained from the Monte Carlo simulation 
and also exact results on a ring of length L=500 for vi = 0.6, a = 0.2 and p = 0.6. The phase transition takes place at 
p c ~ 0.66. 

values of (3. Simulation data also show that this behaviour is generic, at least for the the family of distribution 
functions given by (|34| ) . In figure 6 the speed impurity obtained from the Monte Carlo simulation is plotted as 
a function of f3 for a = 0.2 and a = 0.5. As can be observed the speed of impurity is a decreasing function of 
(3 up to a critical point f3 c . Analytical calculations predict j3 c = 1 — u m — 0.6611834. Above the critical point 
the speed of the impurity is a constant independent of f3. This agrees with our previous results (|2^). As we 
mentioned above, the speed of impurity, as well as the total current of the ordinary particles in the shock phase 
are linear functions of a or f3. In order to see this, we have plotted Jtotai as a function of (3 for a — 0.2 obtained 
from the Monte Carlo simulation and compared with exact result ( p9| ) in figure 7. As can be seen for a fixed 
value of a, Jtotai is a linear function of /3 for (3 < [3 C ~ 0.66 while it remains constant above this point. Now 
for vi < 0.25 we expect that the phase diagram has only two different phases. We take v\ = 0.2 so (3 is limited 
to (3 > 1 — v\ = 0.8. Figure 8 shows the total density profile of the ordinary particles for v\ — 0.2, (3 = 0.85 
and two different values of a. For a fixed value of [3, one can see at most two different phases depending on 
the value of a. In figure 9 we have plotted the speed of the impurity as a function of a for (3 = 0.8, 0.85 and 
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Figure 8: Plot of the total density profile of the ordinary particles on a ring of length L=500 for ,3=0.85, vi—0.2, p = 0.6 
and different values of a. 
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Figure 9: Plot of the speed of the impurity as a function of a obtained from the Monte Carlo simulation on a ring of 
length L=500 for /3 = 0.8, 0.85, 0.9, Vi — 0.2 and, p = 0.6. It is seen that the speed of impurity is independent of j3. 



(3 = 0.9. Analytical results (|29|) predict that for any values of fi the speed of impurity is an increasing function 
of a. As can be seen in figure 9 the speed of impurity is independent of (3 and also a linear function of a for 
a < 0.2. Again this is in agreement with our analytical results (|2^). 
Remark: 

As we mentioned above, the analytical results for the case -4-S > suggest that the phase diagram in the 

u—v\ 

case < a < +oo is made up of two different phases given by (|2^) and (|30|). The Monte Carlo simulations 
for this case have also been carried out. Our results confirm both the phase transition point a c — v\ and also 
the behaviours of V and Jtotai obtained analytically. In figure 10 we have plotted the speed of impurity as a 
function of a for (3 = 0.8, 0.85 and (3 = 0.9. As can be seen the speed of impurity is an increasing function of a 
up to the point a ~ 0.2 which is exactly the transition point in this case. Above the transition point the speed 
of impurity V is no longer a function of a, as the analytical results had predicted. The density profile of the 
ordinary particles on the ring ptotai{i) in the second phase (a > vx) is the same as the phase (|l5|). 
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Figure 10: Plot of the speed of the impurity as a function of a for < a < 1. The results are obtained from the Monte 
Carlo simulation on a ring of length L=500 for f3 — 0.8, 0.85, 0.9, v\ — 0.2 and p = 0.6. It is seen that the speed of 
impurity is independent of f) for a > a c ~ 0.2. 

4 Conclusion 

In this paper we considered the multi-species ASEP (first introduced in p^|) with periodic boundary condition 
in the presence of a single impurity. Using the MPA, we studied the stationary state properties of this model 
by introducing a grand canonical ensemble. We obtained not only the exact phase structure but also exact 
expressions for the mean speed of the impurity V and the total current of the ordinary particles Jtotai in the 
thermodynamic limit. The Monte Carlo Simulation data confirms the correctness of our analytical calculations 



at least for a large family of the density distribution functions given by (34). 

Depending on the sign of (|l2|), which is only a function of the density distribution of the ordinary particles 
on the ring, the phase diagram of the model can exist in two different regimes. Both the analytical and the 
Monte Carlo simulation data show that, as far as the sign of (|lj) is negative, the phase diagram of our model 
has nearly the same structre of the phase diagram of the one-species model Jig. However, there are some 



differences: in comparison to the one-species ASEP the speed of impurity in the shock phase (16), which is in 
fact the speed of the shock in this phase, does not have the simple form a — /3; instead it depends on the density 
distribution function of the ordinary particles on the ring. Also in the multi-species case, the boundaries of 
the phase diagram in the case -^S 



< are specified by the quantity u m which is the solution of (fL9h and 



its value depends on p(v). In the special limit p(v) — v5(v — 1) all our results converge on those obtained in 
fyl . All of the analytical results as well as the Monte Carlo simulation data confirm that if the sign of (|l2|) 
is positive, the phase diagram has no longer four different phases; instead it has one phase. The impurity, by 
definition, should be the slowest particle in the ensemble of the particles on the ring; however, if one releases 
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this condition then a new phase appears. The phase transition point a c — V\ can be calculated analytically and 
also confirmed by simulation. As we mentioned, this new regime cannot be seen in the one-species ASEP. 
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